Real-time polarization mode dispersion characterization

ABSTRACT

In optical media (e.g., an optical fiber) having first and second order polarization mode dispersion, the relationship between the polarization mode dispersion and the change in the measured averaged state of polarization of an optical signal in the media as a function of the signal&#39;s bandwidth can be predictable. Accordingly, first and second order polarization can be determined by measuring a series of the averaged state of polarization&#39;s for an optical signal with varying signal bandwidths.

CROSS-REFERENCE TO RELATED APPLICATIONS

[0001] This application claims priority to U.S. Provisional Application Serial No. 60/327,150, entitled “FIRST AND SECOND ORDER PMD CHARACTERIZATION USING BANDWIDTH-AVERAGED STATE-OF-POLARIZATION” and filed Oct. 4, 2001, which is incorporated herein by reference in its entirety.

BACKGROUND OF THE INVENTION

[0002] The invention relates generally to fiber optics and more specifically to methods of characterizing the polarization mode dispersion in a fiber optic transmission system.

[0003] Polarization mode dispersion (“PMD”) generally refers to variations in the time delay of a polarized optical signal traveling through an optical transmission system, such as a single-mode optical fiber. PMD arises in an optical fiber because of asymmetries in the optical fiber core, such as, for example, core ellipticity created during optical manufacturing and bending stresses resulting from handling or formed during installation. The asymmetries of the fiber core cause random changes in the state of polarization (“SOP”) of optical signals propagating through the fiber. Different SOPs propagate through the optical fiber core at different relative speeds, for example, some SOPs travel faster and some travel slower, resulting in a pulse width distortion of a transmitted optical signal relative to an input optical signal. Additionally, the asymmetries of the core are highly susceptible to environmental fluctuations, such as temperature or movement of the fiber, which occur as fast as milliseconds and result in a time varying pulse width distortion of the transmitted optical signal.

[0004] First order PMD refers to a time delay between two orthogonally polarized principal states of polarization (“PSP”). The PSPs form a convenient basis set to describe and characterize each SOP, and evaluate the effects of PMD in the fiber. Using the PSPs as a basis set, each SOP propagating through an optical fiber is represented by a linear combination of the two orthogonally polarized PSPs. The varying pulse width distortion of the SOP is a function of a varying delay between the PSPs. Theoretically, each PSP experiences a time of flight difference through the optical fiber, commonly known as differential group delay (DGD), resulting in a time delay between the two PSPs at the fiber output. The output SOP is represented by a linear combination of the PSPs that are time delayed with respect to each other. A greater time delay between the PSPs corresponds to a larger relative difference between the input SOP pulse width and the output SOP pulse width.

[0005] Optical fibers have a differential group delay (“DGD”) between the two PSPs on the order of 0.1 ps/(km)^(1/2). In older fiber optic cables, such as the cables used in terrestrial networks, the DGD is on the order of 2.0 ps/(km)^(1/2) and results in time delays of about 50 picoseconds for transmission distances of only several hundred kilometers.

[0006] Second order PMD refers to the frequency dependence of the first order PMD. This is related to chromatic dispersion of signals in optical fibers, which arises due to the frequency dependence of the refractive index of the fiber. Second order PMD causes the dispersion of the PSPs, which further distorts the optical signal propagating in the fiber.

[0007] As the demand for faster optical data transmission increases, such as from gigabits per second to terabits per second, optical pulse width distortion due to PMD may become one of the factors limiting data transmission rate.

SUMMARY

[0008] In optical media (e.g., an optical fiber) having first and second order polarization mode dispersion (“PMD”), there is a predictable relationship between the PMD and the change in the measured averaged state of polarization (“SOP”) of an optical signal in the media as a function of the signal's bandwidth. Accordingly, the first and second order PMD is determined by measuring a series of averaged SOPs for varying signal bandwidths. The signal bandwidth can be varied symmetrically or asymmetrically with respect to a particular frequency.

[0009] In one aspect, a method for characterizing polarization mode dispersion in an optical medium includes measuring an averaged state of polarization vector of an optical signal in the optical medium for each of a plurality of spectra, each spectrum having a frequency bandwidth different from the other spectra of the plurality of spectra. The method further includes determining a characteristic of the polarization mode dispersion from an averaged state of polarization vectors and the frequency bandwidths.

[0010] Embodiments of the method can include one or more of the following features.

[0011] At least a portion of the plurality of spectra have a symmetric bandwidth about a predefined frequency. Alternatively, or additionally, at least a portion of the plurality of spectra have an asymmetric bandwidth about a predefined frequency.

[0012] In some embodiments, the method further includes determining a second state of polarization vector from the measured averaged state of polarization vectors for each of a plurality of spectra, each of the second state of polarization vectors having a symmetric bandwidth about a predefined frequency. Alternatively, or additionally, the method can further includes determining additional state of polarization vector from the measured averaged state of polarization vectors for each of a plurality of spectra, each of the additional state of polarization vectors having an asymmetric bandwidth about a predefined frequency.

[0013] Measuring the averaged state of polarization vector of the optical signal can include deriving a test optical signal from the optical signal, optically filtering the test signal, and measuring averaged state of polarization vectors of the filtered test signal. Optically filtering the test signal can include sequentially directing the test signal through each of a plurality of bandpass filters, and measuring an averaged state of polarization vector for each of the bandpass filters. The frequency spectra of the averaged state of polarization vectors can be related to the bandpass filters.

[0014] Determining the characteristic of the polarization mode dispersion can include measuring a perturbation of the averaged state of polarization vectors from a polarization state at a predefined frequency and a parameter related to the bandwidths of the different spectra. The characteristic of polarization mode dispersion can include a measure of first order polarization mode dispersion and/or a measure of the second order polarization mode dispersion (e.g., a four component vector characterizing the first or second polarization mode dispersion). The characteristic of polarization mode dispersion can be determined in real-time.

[0015] The method can be used to characterize polarization mode dispersion in a number of optical media. For example, the optical medium can be an optical fiber, a birefringent crystal, an optical polymer, or a liquid crystal. Another example is an optical delay line.

[0016] Determining the characteristic of polarization mode dispersion can include determining a trajectory of the averaged state of polarization in the Poincaré sphere as a function of frequency bandwidth.

[0017] In another aspect, an optical system includes a variable optical filter positioned to variably filter the transmission bandwidth of an optical test signal to a polarimeter, and a polarimeter configured to receive the filtered optical test signal and to measure the average state of polarization of the optical test signal for each transmission bandwidth. The system also includes a controller configured to determine first and second order polarization mode dispersion parameters of the optical input signal based on the average state of polarization and the transmission bandwidth of the filtered optical test signal. The optical system can be configured to implement the method for characterizing polarization mode dispersion.

[0018] In yet another aspect, an optical system can include an optical fiber configured to transmit an optical signal between a first location and a second location, a polarimeter configured to receive an optical test signal derived from the optical signal, and to measure an averaged state of polarization for different transmitted bandwidths of the optical input signal, and a controller in communication with the polarimeter and configured to determine first and second polarization mode dispersion parameters of the optical signal based on the averaged state of polarization and transmitted bandwidths of the optical test signal.

[0019] Embodiments of the optical system can include one or more of the following features. The optical system can further include an optical fiber configured to transmit an optical signal from a first location to a second location, wherein the optical test signal is derived from the optical signal. The optical system can also include a compensator configured to compensate the optical signal based on the first and second order polarization mode dispersion parameters to produce a compensated output optical signal. The optical signal can be compensated in a feed-forward mode. The system can include a receiver and/or a transmitter. The receiver can be configured to detect the compensated output optical signal at the second location. The transmitter can be configured to launch the optical signal into the optical fiber at the first location. The transmitter can include a polarization scrambler.

[0020] The optical system can include a tap to derive the optical test signal from the optical signal. The tap can be coupled to the optical fiber near the second location, and can direct the optical test signal into a waveguide that is coupled to the variable optical filter and the polarimeter. First and second order polarization mode dispersion parameters of the optical test signal can be related to first and second order polarization mode dispersion parameters of the optical signal. For example, first and second order polarization mode dispersion parameters of the optical test signal can be the same as first and second order polarization mode dispersion parameters of the optical signal.

[0021] Embodiments of the methods or optical system can include one or more of the following advantages. The methods can characterize first and second order PMD real-time, and can be implemented in feed-forward mode. In some embodiments, the methods can be used to estimate first and second order PMD for setting the starting point for feedback PMD algorithms from which to work. The methods can be implemented without disrupting the operation of the optical telecommunications channels they are monitoring. The methods can be readily implemented in existing telecommunications systems, and apparatus for implementing the methods can be installed at the receiver end of a fiber optic system and can utilize components already in the system. Because the techniques measure the SOP of an optical signal averaged over relatively large bandwidths (e.g., more than 1 GHz, such as 5 GHz or more), measurements can be made with high speed and a high signal-to-noise ratio. The techniques disclosed herein can be implemented to characterize the optical properties of other optical components (e.g., birefringent components).

[0022] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the present specification, including definitions, will control. In addition, the apparatus, methods, and examples are illustrative only and not intended to be limiting.

[0023] Other features and advantages of the invention will be apparent from the following detailed description, and from the claims.

DESCRIPTION OF DRAWINGS

[0024]FIG. 1 is a schematic representation of a fiber optic system including components for characterizing the polarization mode dispersion (“PMD”) in the system;

[0025]FIG. 2 is a schematic representation of a Poincaré sphere;

[0026] FIGS. 3(a)-3(d) are schematic representations of the Poincaré sphere showing the motion of averaged state of polarization (“SOP”) measurements as a function of symmetrically broadening bandwidth for a system with negligible second order PMD;

[0027] FIGS. 4(a) and 4(b) are schematic representations of a section of the Poincaré sphere perpendicular to a principle state of polarization (“PSP”) showing the trajectory of averaged SOP for symmetrically and asymmetrically broadening bandwidth about ω₀, respectively;

[0028]FIG. 5 is a schematic representation of a section of the Poincaré sphere perpendicular to the PSP showing a vector {right arrow over (T)}(ω₀)={right arrow over (M)}₁−{right arrow over (M)}₂;

[0029]FIG. 6(a) is a plot showing the normalized transmission as a function of frequency for a series of optical filters used for measuring average SOPs between frequencies ω_(a) and ω_(b); and

[0030]FIG. 6(b) is a schematic representation of a section of the Poincaré sphere perpendicular to the PSP showing the trajectory of the averaged SOP for SOPs acquired using the optical filters of FIG. 6(a).

[0031] Like reference symbols in the various drawings indicate like elements.

DETAILED DESCRIPTION

[0032] Referring to FIG. 1, a fiber optic system 100 includes a transmitter 110, an optical fiber 120, and a receiver 130. During operation, transmitter 110 launches a polarized optical input signal into one end of optical fiber 120. Transmitter 110 includes a light source, such as a laser diode. Optionally, transmitter 110 can include a polarization scrambler, which randomly varies the polarization state of the optical input signal prior to coupling the signal into optical fiber 120. Optical fiber 120 transmits the optical signal to the opposite end of the fiber where it is detected by receiver 130.

[0033] Fiber optic system 100 also includes a polarization mode dispersion (“PMD”) compensator 140 for reducing the PMD of optical signals transmitted through optical fiber 120. PMD compensator 140 reduces the PMD of optical signals by modifying the polarization state of the optical signal based on the measured PMD of a test signal derived from the optical signal at tap 150. The state of polarization (“SOP”) and PMD of the test signal are related to the SOP and PMD of the optical signal in a predictable way. Thus, measuring the PMD of the test signal provides a measure of the PMD of the optical signal emerging from optical fiber 120. The test signal is passed through a variable optical filter 160, which includes multiple filters of different bandwidth. During operation, variable optical filter 160 sequentially filters the test signal with each of the multiple filters, and a polarimeter 170 measures the averaged SOP of each of the filtered test signals. Polarimeter 170 includes a processor (not shown), which determines one or more characteristics of the PMD of the test signal from the measured average SOPs, and adjusts PMD compensator 140 to reduce the PMD of the optical signal in real-time. PMD compensator 140 works in a feed-forward scheme.

[0034] In general, variable optical filter 160 can be any filter or combination of filters that provide a suitable range of transmission frequencies and transmission bandwidths to determine the PMD of the optical test signal. For example, one can place two filters having overlapping, broadband transmission spectra one after another. Varying the center frequencies of their transmission bands varies the bandwidth of a signal transmitted by the filter pair. Another example is to provide a number of Fabry Perot filters, each having a different transmission bandwidth and to do filtering separately. In general, the range of transmission bandwidths used to filter the test signal depends on the DGD of optical fiber 120. Accordingly, variable optical filter 160 should be selected to provide a sufficient maximum bandwidth and range of bandwidths for the particular system in which it is being implemented.

[0035] To determine the PMD of the test signal, polarimeter 170 measures the averaged SOP of the filtered test signal for several different transmission bands of variable filter 160. Due to the finite response time of detectors used in the polarimeter, the polarimeter intrinsically measures the SOP averaged over the pulse duration. Equivalently, in the frequency domain, the polarimeter averages the SOP over the spectral bandwidth of the transmitted filtered test signal. Accordingly, the measured SOP is an averaged SOP.

[0036] Before describing specific implementations for determining the PMD of fiber 120, it can be useful to describe the mathematical formalism used to determine PMD from the acquired data.

[0037] A frequency-dependent SOP is defined in terms of a Stokes vector:

{right arrow over (S)}(ω)=(S ₀ S ₁ S ₂ S ₃)  (1)

[0038] where

S ₀ =E _(x) ² +E _(y) ²  (2)

S ₁ =E _(x) ² −E _(y) ²  (3)

S ₂=2E _(x) E _(y) cos(φ)  (4)

S ₂=2E _(x) E _(y) sin(φ)  (5)

[0039] and E_(x) and E_(y) are the magnitudes of the x and y component electric field complex amplitudes, respectively. φ is the relative phase between the two.

[0040] Referring to FIG. 2, a convenient and intuitive graphical representation of SOPs is a Poincaré sphere 200. The parameters s₁(ω), s₂(ω) and s₃(ω) are represented by the relationship s_(i)(ω)=(S_(i))/(S₀), where i is 1, 2, or 3, and are used to convert the Stokes parameters to corresponding x, y, and z components in a three dimensional Cartesian coordinate system. In a three dimensional Cartesian coordinate system, Poincaré sphere 200 is defined by a set of points containing all possible SOPs. Points on the surface of the sphere represent purely polarized states, while points inside the sphere represent partially polarized SOPs. As seen in FIG. 2, all linear polarization states are located on an equator 210 of sphere 200, while left and right circular polarizations are located at a north pole 220 and a south pole 230 (opposite north pole 220), respectively. All other points represent elliptical polarizations that lie away from equator 210 and poles 220, 230. Each SOP on sphere 200 is identifiable by its latitude, 2φ, and longitude, 2λ, by using equations:

s ₁(ω)=cos (2φ) cos (2λ)  (6)

s ₂(ω)=cos (2φ) cos (2λ)  (7)

s ₃(ω)=sin (2φ)  (8)

[0041] where any two orthogonal SOPs lie directly opposite each other. For example, a linear vertical polarization 212 is 180 degrees away from a linear horizontal polarization 214 on equator 210.

[0042] Referring to FIG. 1, polarimeter 170 measures a weighted average of the Stokes vector ŝ(ω) across the spectrum of the signal. This averaged SOP, {right arrow over (r)}, is given by $\begin{matrix} {\overset{\rightarrow}{r} = {\int{\frac{\omega}{2\quad \pi}{{f(\omega)}}^{2}{\hat{s}(\omega)}}}} & (9) \end{matrix}$

[0043] where |f(ω)|² is the filtered test signal spectrum normalized to ${\int{\frac{\omega}{2\quad \pi}{{f(\omega)}}^{2}}} = 1.$

[0044] In the limit of a narrowband filtered test signal |f(ω)|²→δ(ω−ω₀), and {right arrow over (r)}→ŝ(ω₀), where, ω₀ is the frequency of a nominally monochromatic pulse.

[0045] The PMD Stokes vector, {right arrow over (Ω)}(ω), describes the trajectory of the test signal SOP, ŝ(ω), on the Poincaré sphere, as a function of frequency, and is related to the test signal SOP by $\begin{matrix} {{\frac{}{\omega}{\hat{s}(\omega)}} = {\overset{\rightarrow}{\Omega} \times {{\hat{s}(\omega)}.}}} & (10) \end{matrix}$

[0046] The magnitude of the PMD Stokes vector is the differential group delay (“DGD”), i.e., |{right arrow over (Ω)}|≡τ_(DGD), and its direction, {circumflex over (Ω)}, is the PSP direction.

[0047] Taking a Taylor expansion of ŝ(ω) around ω₀ up to second order, substituting into Eq. (9) and making use of Eq. (10), one obtains an expression for the measured SOP, {right arrow over (r)}, as a function of the parameters Δω_(pulse) and Δω_(pulse) ²: $\begin{matrix} {\overset{\rightarrow}{r} = {{\hat{s}\left( \omega_{0} \right)} + {\Delta \quad {\omega_{pulse}\left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}} + {\frac{1}{2}\Delta \quad {\omega_{pulse}^{2}\left\lbrack {{\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack} + {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}}} \right\rbrack}}}} & (11) \end{matrix}$

[0048] where ${{\Delta \quad \omega_{pulse}} \equiv {\int{\frac{\omega}{2\quad \pi}{{f(\omega)}}^{2}\left( {\omega - \omega_{o}} \right)}}},{{\Delta \quad \omega_{pulse}^{2}} \equiv {\int{\frac{\omega}{2\quad \pi}{{f(\omega)}}^{2}\left( {\omega - \omega_{o}} \right)^{2}\quad {and}\quad {\overset{\rightarrow}{\Omega}}_{\omega}}} \equiv \frac{\overset{\rightarrow}{\Omega}}{\omega}}$

[0049] is the 2_(nd) order PMD. Eq. (11) expresses {right arrow over (r)} as perturbed from the pure polarization state ŝ(ω₀), by a vector, {right arrow over (M)}, which is defined as $\begin{matrix} {\overset{\rightarrow}{M} = {{\overset{\rightarrow}{r} - {\hat{s}\left( \omega_{0} \right)}} = {{\Delta \quad {\omega_{pulse}\left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}} + {\frac{1}{2}\Delta \quad {\omega_{pulse}^{2}\left\lbrack {{\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack} + {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}}} \right\rbrack}}}}} & (12) \end{matrix}$

[0050] Referring to FIGS. 3(a)-(d), a system with purely first order PMD (i.e., {right arrow over (Ω)}_(ω)=0) such as for a birefringent crystal or a polarization maintaining fiber, is shown. In this case, {right arrow over (M)} is a vector lying in a plane normal to the PSP axis. For a monochromatic signal having Stokes vector ŝ(ω₀), {right arrow over (M)}=0 and {right arrow over (r)} is a point on the surface of the Poincaré sphere. If {right arrow over (r)} is measured over a finite bandwidth in a system having first order PMD (i.e., non-zero {right arrow over (Ω)}) as is usually the case in real systems, it represents an averaged SOP that is partially polarized (i.e., |{right arrow over (r)}|<1) and {right arrow over (r)} departs the surface of the Poincaré sphere. As one increases the spectral bandwidth of the signal, the averaged SOP prescribes a trajectory originating on the surface of the Poincaré sphere and moving inward toward the PSP.

[0051] Referring to FIGS. 4(a) and 4(b), for a system with purely first order PMD and a spectrum that is symmetrical about ω₀ (i.e., Δω_(pulse)=0), {right arrow over (M)} reduces to $\begin{matrix} {\overset{\rightarrow}{M} = {\frac{1}{2}\Delta \quad \omega_{pulse}^{2}\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}} & (13) \end{matrix}$

[0052] The {right arrow over (Ω)}×{right arrow over (Ω)}× term is a reduction of components orthogonal to {right arrow over (Ω)} (generally called {right arrow over (s)}⊥(ω₀), given by ${\overset{\rightarrow}{s}\bot} = {\hat{s} - {\frac{\left( {\hat{s} \cdot \overset{\rightarrow}{\Omega}} \right)\overset{\rightarrow}{\Omega}}{\tau_{DGD}^{2}}.}}$

[0053] The following equations result: $\begin{matrix} {\overset{\rightarrow}{r} = {{{\hat{s}\left( \omega_{0} \right)} - {\frac{1}{2}\Delta \quad \omega_{pulse}^{2}\tau_{DGD}^{2}\overset{\rightarrow}{s}}}\bot\left( \omega_{0} \right)}} & (14) \end{matrix}$

[0054] and $\begin{matrix} {\overset{\rightarrow}{M} = {{{- \frac{1}{2}}\Delta \quad \omega_{pulse}^{2}\tau_{DGD}^{2}\overset{\rightarrow}{s}}\bot{\left( \omega_{0} \right).}}} & (15) \end{matrix}$

[0055] Recalling that {circumflex over (Ω)} is parallel to the PSP, the vector {right arrow over (M)} points toward the center of an arc that is traced out by ŝ(ω) when ω is varied, as seen in FIG. 4(a). Hence, the averaged SOP, {right arrow over (r)}, aligns toward the PSP direction as the spectral bandwidth increases.

[0056] The degree of polarization (“DOP”) of the measured SOP as a function of bandwidth, Δω_(pulse) ², is given by

DOP ² ≡|{right arrow over (r)}| ²=1−Δω_(pulse) ²τ_(DGD) ² {right arrow over (s)}⊥(ω₀)|²  (16)

[0057] Accordingly, Eqs. (14) and (16) provide three independent scalar equations to solve for {right arrow over (Ω)}. The straight-line plots of the components of {right arrow over (r)} and DOP² are linear in Δω_(pulse) ², with gradients ${\overset{\rightarrow}{g}}_{r} = {{{- \frac{1}{2}}\tau_{DGD}^{2}\overset{\rightarrow}{s}}\bot\left( \omega_{0} \right)}$

[0058] and g_(Dop)=−τ_(DGD) ²|{right arrow over (s)}⊥(ω₂)|², respectively. These equations can be rearranged as $\tau_{DGD} = \frac{2\quad {{\overset{\rightarrow}{g}}_{r}}}{\sqrt{g_{DOP}}}$

[0059] and ${\overset{\rightarrow}{s}\bot\left( \omega_{0} \right)} = {\frac{{- 2}{\overset{\rightarrow}{g}}_{r}}{\tau_{DGD}^{2}}.}$

[0060] By extracting these gradients from a series of measurements, τ_(DGD) and {right arrow over (s)}⊥(ω₀) can be determined. As a result of knowing {right arrow over (s)}⊥(ω₀), the PSP direction is known since the unit vector of {right arrow over (Ω)} lies in the same direction of ŝ(ω₀)−{right arrow over (s)}⊥(ω₀). ŝ(ω₀) is obtained from the y-intercepts of Eq. (14).

[0061] The formalism presented above treats the situation where the spectrum is symmetric about ω₀. However, if the spectrum is asymmetric about ω₀, {right arrow over (M)} becomes $\begin{matrix} {\overset{\rightarrow}{M} = {{\Delta \quad {\omega_{pulse}\left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}} + {\frac{1}{2}\Delta \quad \omega_{pulse}^{2}\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}}} & (17) \end{matrix}$

[0062] An increase of bandwidth still moves the measured SOP inward as shown in FIG. 4(b). While the trajectory of the averaged SOP no longer moves toward the PSP axis, it still lies in the plane whose normal vector is in the PSP direction. Thus, {right arrow over (M)} still lies in this plane.

[0063] Second order PMD can also be present. In such systems, {right arrow over (M)} no longer lies in the plane normal to the PSP due to an additional vector component, ${\frac{1}{2}\Delta \quad {\omega_{pulse}^{2}\quad\left\lbrack {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}},$

[0064] in Eq. (12).

[0065] For the case of a symmetric spectrum (i.e. Δω_(pulse)=0), $\begin{matrix} {\overset{\rightarrow}{r} = {{\hat{s}\left( \omega_{0} \right)} + {\frac{1}{2}\Delta \quad {\omega_{pulse}^{2}\left\lbrack {{\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack} + {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}}} \right\rbrack}}}} & (18) \\ {\overset{\rightarrow}{M} = {\frac{1}{2}\Delta \quad {\omega_{pulse}^{2}\left\lbrack {{\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack} + {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}}} \right\rbrack}}} & (19) \end{matrix}$

[0066] From Eqs. (18) and (19), it is apparent that taking the averaged SOP over a symmetrical spectrum, the 2^(nd) order PMD contributes to the same order of magnitude as 1^(st) order PMD in {right arrow over (r)} and {right arrow over (M)}. An estimate of the 2^(nd) order PMD parameters can be obtained using the following methodology: (i) isolate the 1^(st) order PMD effect, determine {right arrow over (Ω)}, and (ii) determine the 2^(nd) order PMD, {right arrow over (Ω)}_(ω), using Eq. (19).

[0067] To isolate 1^(st) order PMD effects, two different signal spectra, |f₁|² and |f₂|², can be used, which give different bandwidths {Δω₁, Δω₁ ²} and {Δω₂, Δω₂ ²} respectively. The corresponding {right arrow over (M)} vectors, from Eq. (12), are: $\begin{matrix} {{\overset{\rightarrow}{M}}_{1} = {{\Delta \quad {\omega_{1}\left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}} + {\frac{1}{2}\Delta \quad {\omega_{1}^{2}\left\lbrack {{\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack} + {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}}} \right\rbrack}}}} & (20) \end{matrix}$

$\begin{matrix} {{\overset{\rightarrow}{M}}_{2} = {{\Delta \quad {\omega_{2}\left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack}} + {\frac{1}{2}\Delta \quad {\omega_{2}^{2}\left\lbrack {{\overset{\rightarrow}{\Omega} \times \left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack} + {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{0} \right)}}} \right\rbrack}}}} & (21) \end{matrix}$

[0068] By multiplying Eq. (20) by Δω₂ ², Eq. (21) by Δω₁ ², and subtracting one from the other, one removes the second order PMD dependence. Defining, ${{\overset{\rightarrow}{T}\left( \omega_{0} \right)} = {{\frac{\Delta \quad \omega_{2}^{2}}{\Delta \quad \omega_{1}^{2}}{\overset{\rightarrow}{M}}_{1}} - {\overset{\rightarrow}{M}}_{2}}},$

[0069] a vector is obtained that depends only on 1st order PMD: $\begin{matrix} {{\overset{\rightarrow}{T}\left( \omega_{0} \right)} = {\left\lbrack {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{0} \right)}} \right\rbrack \left( {{\Delta \quad \omega_{1}\frac{\Delta \quad \omega_{2}^{2}}{\Delta \quad \omega_{1}^{2}}} - {\Delta \quad \omega_{2}}} \right)}} & (22) \end{matrix}$

[0070] Using different sets of spectra {|f₁|²,|f₂|²}, one deduces {right arrow over (Ω)}×ŝ(ω₀) from the gradients of Eq. (22) when various components of vector {right arrow over (T)}(ω₀) are plotted with various sets of $\left( {{{\Delta\omega}_{1}\frac{{\Delta\omega}_{2}^{2}}{{\Delta\omega}_{1}^{2}}} - {\Delta\omega}_{2}} \right).$

[0071] It is worthwhile to note that {right arrow over (T)}(ω₀) is parallel {right arrow over (Ω)}×ŝ(ω₀), which is the tangential vector to the arc on the Poincaré sphere traced out by ŝ(ω₀) when ω is varied. Recalling that ${{\frac{\quad}{\omega}{\hat{s}(\omega)}} = {\overset{\rightarrow}{\Omega} \times {\hat{s}(\omega)}}},$

[0072] alternatively {right arrow over (Ω)}×ŝ(ω₀) can be approximated from the differential change of the Stokes vector, ŝ(ω), when ω is changed by an infinitesimal amount. However, the aforementioned method is more accurate since one estimates {right arrow over (Ω)}×ŝ(ω₀) from the gradients of straight-line plots of Eq. (22) with various bandwidths of finite-bandwidth spectra, and, by virtue that averaged SOP is measured.

[0073] Referring to FIG. 5, to help visualize {right arrow over (T)}(ω₀) and its relationship with {right arrow over (M)}₁ and {right arrow over (M)}₂, a special case in which a combination of asymmetrical |f₁|² and symmetrical |f₂|², where Δω₂ ²=Δω₁ ² are used. Because |f₂|² is symmetrical, Δω₂=0. Thus, {right arrow over (T)}(ω₀) reduces to

{right arrow over (T)}(ω₀)={right arrow over (M)} ₁ −{right arrow over (M)} ₂=Δω₁ {right arrow over (Ω)}×ŝ(ω₀)  (23)

[0074] To get {right arrow over (Ω)} from {right arrow over (Ω)}×ŝ(ω₀), knowledge of another vector, for example, {right arrow over (Ω)}×ŝ′, can be needed. By taking the cross-product of {right arrow over (Ω)}×ŝ(ω₀) and {right arrow over (Ω)}×ŝ′, and using the vector identity a×(b×c)=(a.c)b−(a.b)c, the first order PMD is obtained from the relationship $\begin{matrix} {\overset{\rightarrow}{\Omega} = \frac{\left( {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{o} \right)}} \right) \times \left( {\overset{\rightarrow}{\Omega} \times \hat{s}} \right)}{\left( {\overset{\rightarrow}{\Omega} \times {\hat{s}\left( \omega_{o} \right)}} \right) \cdot \hat{s}}} & (24) \end{matrix}$

[0075] There are at least two different ways to experimentally obtain {right arrow over (Ω)}×ŝ′:

[0076] The first method makes use of the bandwidth of the PSP. For an optical frequency ω₀ lying within this PSP bandwidth around ω₀, one can approximate {right arrow over (Ω)}(ω₁)≈{right arrow over (Ω)}(ω₀)={right arrow over (Ω)}. So, by repeating the above measurement and analysis and using Eq. (22) about ω₁, {right arrow over (Ω)}×ŝ(ω₁) is obtained for the calculation of {right arrow over (Ω)}. It is worthwhile to note if ω₀ and ω₁ are separated but stay within this PSP bandwidth, the angle subtended by {right arrow over (Ω)}×ŝ(ω₀) and {right arrow over (Ω)}×ŝ(ω₁) is quoted as 0.7 radian by C. D. Poole et al., in “Polarization Dispersion and Principle States in a 147 km Undersea Lightwave Cable,” J. Lightwave Tech., Vol. 6, No. 7, pp 1185-1190 (July, 1988), and as 0.8±0.4 radian by S. Betti et al., in “Evolution of the bandwidth of the principle states of polarization in single mode fibers,” Opt. Lett., Vol. 16, No. 7, pp 467-469 (April, 1991), respectively. Both these references are hereby incorporated by reference in their entirety.

[0077] The second method uses another set of measurements made using a second input SOP.

[0078] By repeating the measurement and analysis for the second input SOP, {right arrow over (Ω)}×ŝ′ is obtained. This scheme can be feasible in optical communication systems where a low speed polarization scrambler (e.g., ˜kHz) is implemented, e.g., to overcome polarization hole burning in an erbium-doped fiber amplifier (EDFA). Polarization scrambling provides multiple different input SOPs for PMD characterization every millisecond interval, which can improve characterization accuracy. By using a least square fit algorithm to find the unit vector that is perpendicular to all these {right arrow over (Ω)}×ŝ′ vectors, one can determine the orientation of the PSP.

[0079] Second order PMD, {right arrow over (Ω)}_(ω), is determined according to Eq. (19). The information listed in Eq. (19) can be determined from a plot of components of {right arrow over (M)} vs. Δω_(pulse) ², which is substantially linear for data acquired using the symmetrical spectrum, |f₂|². The empirically determined gradients, {right arrow over (g)}, of these plots correspond to the components of {right arrow over (Ω)}×[{right arrow over (Ω)}×ŝ(ω₀)]+{right arrow over (Ω)}_(ω)×ŝ(ω₀). Hence, the components of {right arrow over (Ω)}_(ω)×ŝ(ω₀) are determined by substituting the known components of {right arrow over (Ω)} and ŝ(ω₀) into the right hand side of

{right arrow over (Ω)}_(ω) ×ŝ(ω ₀)={right arrow over (g)}−{right arrow over (Ω)}×[{right arrow over (Ω)}×ŝ(ω₀)]  (25)

[0080] Using a method similar to the aforementioned methods for obtaining {right arrow over (Ω)}×ŝ′, the vector {right arrow over (Ω)}_(ω) ×ŝ′. By taking the cross-product of {right arrow over (Ω)} _(ω) ×ŝ(ω ₀) and {right arrow over (Ω)}_(ω) ×ŝ′, applying the vector identity a×(b×c)=(a.c)b−(a.b)c, the second order PMD can be determined from the equation: $\begin{matrix} {{\overset{\rightarrow}{\Omega}}_{\omega} = \frac{\left( {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{o} \right)}} \right) \times \left( {{\overset{\rightarrow}{\Omega}}_{\omega} \times \hat{s}} \right)}{\left( {{\overset{\rightarrow}{\Omega}}_{\omega} \times {\hat{s}\left( \omega_{o} \right)}} \right) \cdot \hat{s}}} & (26) \end{matrix}$

[0081] There are a number of ways to measure the SOP of a signal for symmetric and asymmetric bandwidth broadening and to implement the aforementioned methods. Referring to FIGS. 6(a) and 6(b), one scheme can include measuring averaged SOPs using several different optical filters that have transmission curves as shown in FIG. 6(a). The unnormalized transmitted spectra (|F₁|² . . . |F_(N)|²) are the products of the on-line signal spectra and these transmission curves. Optical frequencies, ω_(a) and ω_(b), in FIG. 6(a) can be separated by as much as about twice the PSP bandwidth. FIG. 6(b) shows the corresponding trajectory 610 of the averaged SOP when the transmission bandwidth of the filters is increased.

[0082] From these measured averaged SOPs, the averaged SOP for the difference of two overlapping spectra can be computed. For example, if the average SOP {right arrow over (r)}₁ for spectrum |F₁|² and {right arrow over (r)}₂ for spectrum |F₂|² are known, the averaged SOP for the difference between these spectra, (|F₂|²−|F₁|²) can be computed using the equation $\begin{matrix} {{\overset{\rightarrow}{r}}_{diff} = \frac{{{\overset{\rightarrow}{r}}_{2}\left( {\int\left| F_{2} \middle| {}_{2}{\omega} \right.} \right)} - {{\overset{\rightarrow}{r}}_{1}\left( {\int\left| F_{1} \middle| {}_{2}{\omega} \right.} \right)}}{\int\left| F_{2} \middle| {}_{2}{{\omega} - \int} \middle| F_{1} \middle| {}_{2}{\omega} \right.}} & (27) \end{matrix}$

[0083] Similarly, the averaged SOP of the sum of two non-overlapping spectra can be computed.

[0084] For example, if one knows {right arrow over (r)}_(diff) for (|F₂|²−|F₁|²) and {right arrow over (r)}′_(diff) for (|F₃|²−|F₂|²), one can obtain the averaged SOP for the combined spectra of (|F₂|²−|F₁|²) and (|F₃|²−|F₂|²) spectra by $\begin{matrix} {{\overset{\rightarrow}{r}}_{sum} = \frac{\begin{matrix} {{{\overset{\rightarrow}{r}}_{diff}\left( {\int{\left( \left| F_{2} \middle| {}_{2}{- \left| F_{1} \right|^{2}} \right. \right){\omega}}} \right)} +} \\ {{\overset{\rightarrow}{r}}_{diff}\left( {\int{\left( \left| F_{3} \middle| {}_{2}{- \left| F_{2} \right|^{2}} \right. \right){\omega}}} \right)} \end{matrix}}{\left( {\int{\left( \left| F_{2} \middle| {}_{2}{- \left| F_{1} \right|^{2}} \right. \right){\omega}}} \right) + \left( {\int{\left( \left| F_{3} \middle| {}_{2}{- \left| F_{2} \right|^{2}} \right. \right){\omega}}} \right)}} & (28) \end{matrix}$

[0085] Using Eqs. (27) and (28) it is possible to determine SOPs corresponding to the sum or difference of any combination of measured spectra.

[0086] To characterize the first and second order PMD of a fiber using the methodology outlined above, four sets of averaged SOP data can be determined: one for symmetrical bandwidth broadening and one for asymmetrical bandwidth broadening of bandwidths about two different frequencies, ω₀ and ω₁, where ω_(a)<ω₀<ω₁<ω_(b). These four data sets are constructed by subtracting and combining spectra of (|F₁|² . . . |F_(N)|²) using Eqs. (27) and (28). Once constructed, the processor initially determines the first order PMD using Eq. (24), and then determines the second order PMD using Eq. (26).

[0087] It will be understood that various modifications may be made to the present embodiment without departing from the spirit and scope of the method and system. For example, although the described system is used to characterize PMD in an optical fiber, the techniques disclosed herein can be used to characterize other optical media. For example, these techniques can be used to determine retardation (i.e., first order PMD) and/or birefringence of a birefringent crystal, polymer, or liquid crystal. Another example includes characterizing the DGD of an optical delay line. Accordingly, other embodiments are within the scope of the following claims. 

What is claimed is:
 1. A method for characterizing polarization mode dispersion in an optical medium, the method comprising: measuring an averaged state of polarization vector of an optical signal in the optical medium for each of a plurality of spectra, each spectrum having a frequency bandwidth different from the other spectra of the plurality of spectra; and determining a characteristic of the polarization mode dispersion from an averaged state of polarization vectors and the frequency bandwidths.
 2. The method of claim 1, wherein at least a portion of the plurality of spectra have a symmetric bandwidth about a predefined frequency.
 3. The method of claim 1, wherein at least a portion of the plurality of spectra have an asymmetric bandwidth about a predefined frequency.
 4. The method of claim 1, further comprising determining a second state of polarization vector from the measured averaged state of polarization vectors for each of a plurality of spectra, each of the second state of polarization vectors having a symmetric bandwidth about a predefined frequency.
 5. The method of claim 1, further comprising determining a second state of polarization vector from the measured averaged state of polarization vectors for each of a plurality of spectra, each of the second state of polarization vectors having an asymmetric bandwidth about a predefined frequency.
 6. The method of claim 1, wherein measuring the averaged state of polarization vector of the optical signal comprises deriving a test optical signal from the optical signal, optically filtering the test signal, and measuring averaged state of polarization vectors of the filtered test signal.
 7. The method of claim 6, wherein optically filtering the test signal comprises sequentially directing the test signal through each of a plurality of bandpass filters, and measuring an averaged state of polarization vector for each of the bandpass filters.
 8. The method of claim 1, wherein determining the characteristic of the polarization mode dispersion includes measuring a perturbation of the averaged state of polarization vectors from a polarization state at a predefined frequency and a parameter related to the bandwidths of the different spectra.
 9. The method of claim 1, wherein the characteristic of polarization mode dispersion comprises a measure of first order polarization mode dispersion.
 10. The method of claim 1, wherein the characteristic of polarization mode dispersion comprises a measure of second order polarization mode dispersion.
 11. The method of claim 1, wherein the characteristic of polarization mode dispersion is determined in real-time.
 12. The method of claim 1, wherein the optical medium is an optical fiber.
 13. The method of claim 1, wherein determining the characteristic of polarization mode dispersion comprises determining a trajectory of the averaged state of polarization in the Poincaré sphere as a function of frequency bandwidth.
 14. An optical system comprising: a variable optical filter positioned to variably filter the transmission bandwidth of an optical test signal to a polarimeter; and a polarimeter configured to receive the filtered optical test signal and to measure the average state of polarization of the optical test signal for each transmission bandwidth; and a controller configured to determine first and second order polarization mode dispersion parameters of the optical input signal based on the average state of polarization and the transmission bandwidth of the filtered optical test signal.
 15. The optical system of claim 14, further comprising an optical fiber configured to transmit an optical signal from a first location to a second location, wherein the optical test signal is derived from the optical signal.
 16. The optical system of claim 14, further comprising a compensator configured to compensate the optical signal based on the first and second order polarization mode dispersion parameters to produce a compensated output optical signal.
 17. The optical signal of claim 16, wherein the optical signal is compensated in a feed-forward mode.
 18. The optical system of claim 16, further comprising a receiver configured to detect the compensated output optical signal at the second location.
 19. The optical system of claim 14, further comprising a transmitter configured to launch the optical signal into the optical fiber at the first location.
 20. The optical system of claim 19, wherein the transmitter comprises a polarization scrambler.
 21. The optical system of claim 14, further comprising a tap configured to derive the optical test signal from the optical signal.
 22. The optical system of claim 21, wherein the tap is coupled to the optical fiber near the second location.
 23. The optical system of claim 14, wherein first and second order polarization mode dispersion parameters of the optical test signal are related to first and second order polarization mode dispersion parameters of the optical signal.
 24. An optical system, comprising: an optical fiber configured to transmit an optical signal between a first location and a second location; a polarimeter configured to receive an optical test signal derived from the optical signal, and to measure an averaged state of polarization for different transmitted bandwidths of the optical input signal; and a controller in communication with the polarimeter and configured to determine first and second polarization mode dispersion parameters of the optical signal based on the averaged state of polarization and transmitted bandwidths of the optical test signal. 